Break-up fragment topology in statistical multifragmentation models 
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Break-up fragmentation patterns together with kinetic and configurational energy fluctuations 
are investigated in the framework of a microcanonical model with fragment degrees of freedom over 
a broad excitation energy range. As far as fragment partitioning is approximately preserved, energy 
fluctuations are found to be rather insensitive to both the way in which the freeze-out volume is 
constrained and the trajectory followed by the system in the excitation energy - freeze-out volume 
space. Due to hard-core repulsion, the freeze-out volume is found to be populated un-uniformly, its 
highly depleted core giving the source a bubble-like structure. The most probable localization of 
the largest fragments in the freeze-out volume may be inferred experimentally from their kinematic 
properties, largely dictated by Coulomb repulsion. 

PACS numbers: 25.70.Mn, 25.70.Pq, 24.10.Pa, 64.60.an 



I. INTRODUCTION 

For more than two decades, nuclear multifragmenta- 
tion benefits from a constant scientific interest whose 
main motivation is the observation of a (liquid-gas-like) 
phase transition at sub-atomic scale 0, . 

Relying on the presumptive existence of an equili- 
brated break-up stage in the simultaneous multi-particle 
decay of excited nuclei, statistical models with cluster de- 
grees of freedom @, 0, IE H 0| represent particularly useful 
tools for the characterization of the equilibrated state of 
the source and, not less important, the study of the as- 
sociated thermodynamics. The remarkable advantage of 
realistically incorporating most properties of bound and 
continuum states via empirical paramctcrizations of clus- 
ter energies or level densities explains their ability to well 
describe a wealth of experimental data produced over a 
broad energetic domain. 

It was demonstrated that experimental data corre- 
sponding to a well-defined equilibrated source may be 
described by a unique solution of such a statistical model 
[1, Q . It is nevertheless not true that the different sta- 
tistical models converge to the same equilibrated source 
if the analysis is done by exclusively considering experi- 
mental (after-burner) information [1 01 ] . This is partly due 
to the different thermodynamical constraints imposed to 
the employed statistical ensembles or mathematical tricks 
designed in order to simplify the partition function or 
speed-up the simulation and, to a much larger extend, to 
the differences in the break-up fragment definition. 

The aim of the present work is to contribute to a deeper 
understanding of the break-up stage of the multifragmen- 
tation decay as ruled by statistical laws. For this reason, 
contributions from dynamics (as radial collective flow) 
and sequential particle evaporations from primary frag- 
ments will be referred to only tangentially, despite that 
over an important region of the considered energy domain 
they play an important role. For the same reason we will 
ignore also eventual fragments recombination subsequent 



to the break-up, thoroughly considered by some authors 
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• whether fluctuations of different energetic degrees 
of freedom are mainly dictated by the localization 
of the decay event into the phase diagram or, con- 
versely, by the dominant fragmentation modes, 

• whether break-up nuclear matter distribution is 
uniform and, if not, 

• whether is it possible to trace the un-homogencities 
from experimentally accessible data. 

The paper is organized as follows: Sec. II offers a 
brief review on the statistical models of multifragmen- 
tation with a special focus on the microcanonical ones, 
employed here; Sec. Ill investigates the sharing of sys- 
tem's available energy among different degrees of freedom 
and the sensitivity of the energy fluctuations to the sys- 
tem phase properties and fragment partition; Sec. IV 
focuses on break-up patterns and the extend in which 
these may be inferred from kinetic energy distributions. 
Modifications of fragment charge distributions brought 
by considering that, at variance with the standard break- 
up picture, primary fragments interact through nuclear 
forces are also addressed in Section IV. Conclusions are 
drawn in Sec. V. 



II. STATISTICAL TREATMENT OF 
MULTIFRAGMENTATION 

Under the equilibrium hypothesis, statistical models 
reduce the physical problem under study to the estima- 
tion of the number of microscopic states compatible with 
the thermodynamical macroscopic constraints. This im- 
plies that assuming that it is possible to write down the 
mathematical expression of the statistical weight of a 
configuration Wc in the appropriate statistical ensem- 
ble, all the thermodynamic quantities may be calculated 
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out of the characteristic partition sum, 



(1) 



while any ensemble-averaged observable X may be ex- 
pressed as, 



<X C >-- 



Ec w cX c 
Ec w c 



(2) 



While for relatively large extensive systems, thermody- 
namical properties are not sensitive to the way in which 
the statistical ensemble is defined, when dealing with 
small systems, as the nuclear ones, it is important to 
choose the most appropriate replica of the physical phe- 
nomenon. The lack of any thermal or chemical potential 
reservoirs in the case of isolated multifragmenting nuclei, 
recommend the microcanonical ensemble as the most rea- 
sonable choice 

HUSH- In this 

case, it is obvious that 
the conserved quantities are the total protons (Z) and 
neutrons (A — Z) numbers, the total energy (E), total 
momentum (P) and, eventually, total angular momen- 
tum (L). The freeze-out volume (V) may be considered 
either as fixed, either as fluctuating. 

Defining a generic break-up configuration by the iso- 
topic, internal and translational properties of each frag- 
ment, C = {A 1 ,Z 1 ,ei,r 1 ,...,A Nc ,Z Nc ,e Nc ,r Nc }, one 
gets for the statistical weight of the constant volume en- 
semble the equation Q, 



W C (A,Z,E,V) oc ^ttfj (^( mAn )^ 
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where I is the moment of inertia, K is the thermal ki- 
netic energy and f2 = xV N ° stays for the free volume or, 
cquivalently, accounts for inter-fragment interaction in 
the hard-core idealization. From Eq. ([3]) it is straightfor- 
ward to calculate the statistical weight of a microcanon- 
ical ensemble with fluctuating volume as, 



W C (A,Z,E,\) 



J W c {A,Z,E,V)exp(-XV)dV. (4) 



It is worthwhile to mention at this point that working 
under a fixed total energy constraint, it results that the 
thermal kinetic energy, a key thermodynamic quantity 
related to the temperature through T -1 = (dS/dE) = 
l/W(A,Z,E,V)dW(A,Z,E,V)/dE, is determined by 
the amount of the energy available after extracting from 
the source excitation the costs of fragment formation 
X^j-Bj, fragment internal excitation J^. e, and mutual 
fragment interaction (X^<j ^j')> 



K = E P 



q-E^-E^ 



(5) 
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FIG. 1: (Color online) Heat capacity versus source excita- 
tion energy for the multifragmenting nucleus (130,60) which 
evolves through the phase space following (1) a constant 
freeze-out volume V — 6Vb path, or, (2) a path character- 
ized by an average freeze-out volume which increases with 
the excitation as indicated in the inset. 



This implies that also the fluctuations of K are strongly 
dependent on the fluctuations of the other three energetic 
degrees of freedom, as we shall see later on. 

The results discussed hereafter have been obtained in 
the framework of the Microcanonical Model of Multifrag- 
mentation (MMM) Q in the case of the medium size nu- 
cleus (130,60) within the commonly accepted scenario ac- 
cording to which the break-up fragments do not interact 
otherwise than via Coulomb forces. The consequences of 
considering in the spirit of Rcfs. 
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up fragments feel also the nuclear proximity potential arc 
discussed only with respect to fragment charge distribu- 
tions, for the sake of completeness. Despite the particular 
choices regarding the model and the nucleus, the results 
are considered generic for the statistical break-up of mul- 
tifragmenting nuclei. 

Two arbitrary paths in the phase diagram have been 
considered, a constant volume path V = 6Vq (full sym- 
bols) and one along which the average volume increases 
with excitation energy (open symbols). The motivation 
of choosing a constant volume path is twofold. On one 
hand, it reproduces the fixed freeze-out volume statistical 
constrained which was used for treating multifragmenta- 
tion over almost two decades and, on the other hand, it 
accounts for the belief that the freeze-out volume (aver- 
age) value docs not change significantly while increasing 
source excitation energy. Twofold is also the motivation 
of choosing the second path. First, it cancels the statis- 
tical constraint of constant volume and, secondly, it ac- 
counts for a freeze-out volume whose (average) value may 
increase with energy, as recent analyses of experimen- 
tal data indicate 



17| . In this last case, the average 
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freeze-out volume increases from 3.5Vb at 2 MeV/nucleon 
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to about 10.4Vb at 14 MeV/nucleon, as indicated in the 
inset of Fig. [1] Even more importantly, the two paths 
differ by the regions of the system phase diagram they 
explore. Thus, following the evolution of the heat capac- 
ity, 



cr 
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(6) 



plotted in Fig. [T] as a function of excitation energy, 
one may notice that the constant volume path is supra- 
critical, while the increasing average-volume path crosses 
the phase coexistence region. Phase coexistence is sig- 
naled by negative values of the heat capacity. 



III. ENERGY SHARING AT BREAK-UP: 
AVERAGE AND RMS VALUES 

Fig. [H presents the average values of total frag- 
ment binding energy (upper panel), internal excitation 
(second upper panel), Coulomb interaction (third up- 
per panel) and thermal kinetic energy (bottom panel) 
(left column) together with their RMS values (right col- 
umn) corresponding to the break-up stage of a (130,60) 
nucleus whose excitation energy ranges from 2 to 14 
MeV/nucleon along the two considered trajectories. Dis- 
tributions of the mean charge of the largest fragment and 
its RMS arc superimposed on the third upper panels with 
full and open stars. Dashed lines are used in the bottom 
panels to indicate how total fragment kinetic energy, a 
quantity experimentally accessible, behaves with respect 
to source excitation. 

One can see that, irrespectively the considered path, 
the more and more advanced fragmentation allowed by an 
increasing source energy, suggested by a rapidly decreas- 
ing Z max , leads to a monotonic diminish of the total bind- 
ing energy and a monotonic increase of the total Coulomb 
interaction energy. The total binding energy decrease is 
due to the increasing fragment surfaces while the increase 
of the total Coulomb interaction energy is explained by 
an increasingly uniform occupation of the volume. While 
the curves corresponding to the two considered paths di- 
verge with excitation, they are still not far one from an- 
other, as for a given E ex their values differ by at most 
20% in the considered energy domain. In contrast with 
this, the amount of energy dissipated in fragment internal 
excitation has a more complex evolution and the relative 
difference among the values obtained along the two paths 
reaches 50% at E ex =14 MeV/nucleon. Nevertheless, the 
evolution and relative magnitude of the above quantities 
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FIG. 2: (Color online) Evolution with source excitation en- 
ergy of mean (left column) and RMS (right column) values 
of total binding energy (upper panel) , total fragment internal 
excitation (second upper panel), Coulomb fragment-fragment 
interaction (third panel) and thermal kinetic energy (bottom 
panel) corresponding to the break-up stage of the (130, 60) 
multifragmenting nucleus. The considered states along a con- 
stant volume, V — 6Vb, (1) and an average volume increasing 
with excitation (2) paths are represented with solid and, re- 
spectively, open circles. Third raw panels depict also, with 
open and solid stars, the excitation energy dependence of the 
largest fragment in each decay event. The dashed lines in the 
bottom panel stand for total fragment kinetic energy. 



are such that the kinetic energy increases monotonically, 
as one would expect (see left bottom panel). 

The right panels of Fig. [2] present the energy fluc- 
tuations and indicate that, as more and more fragment 
partitions are possible with the increasing energy, a(B) 
and a{E in t) rise as well. Very interestingly, aiVc) aug- 
ments up to 4 MeV/nucleon and then decreases. The 
positive slope region corresponds to the energy domain 
where configurations containing one heavy residue are 
dominant. The negative slope interval corresponds to a 
regime of rather advanced fragmentation which allows for 
a more uniform population of the freeze-out volume. As 
one may notice, the peak of cr(Vc) corresponds roughly 
to the peak of cr(Z max ) (full and open stars in the third 
right panel) and indicates that the largest fragment Z max 
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dictates the geometrical arrangement of fragments and, 
finally, the Coulomb energy. 

Another remark is that, because of the fact that the re- 
duction of Coulomb energy fluctuation is less significant 
than the increase of internal excitation and binding en- 
ergy fluctuations, a(K t h) increases monotonically. Nev- 
ertheless, analysing the experimentally accessible frag- 
ment total kinetic energy distribution, one would note a 
peak at 3 MeV/nucleon, as the consequence of summing 
up the peaked u(Vc) with the monotonically increasing 
<j(K th ) (right bottom panel). 

But, the first important result is that fluctuations of ki- 
netic and configurational energetic channels prove rather 
insensitive to both freeze-out volume constraints and the 
trajectory followed by the system into the excitation en- 
ergy - freeze-out volume plane, provided that the frag- 
mentation pattern is preserved. The result is striking the 
more as the two considered trajectories explore different 
regions of the phase diagram. 



IV. FRAGMENTATION PATTERNS AND 
NUCLEAR MATTER RADIAL DISTRIBUTIONS 

The break-up fragmentation pattern corresponding to 
4 MeV/nucleon excitation energy, where the largest fluc- 
tuations in cr(Z max ) and cr(Vc) manifest themselves, is 
illustrated in the upper panel of Fig. [3] while the up- 
per panel of Fig. [4] presents the fragmentation pat- 
tern obtained at a slightly higher source excitation, 6 
MeV/nucleon. As no sensitivity was found to the way in 
which the frcczc-out volume in constrained, from here on 
we shall consider only the case corresponding to V = 6Vq. 

One can see that at E ex = 4 MeV/nucleon the dom- 
inant fragmentation mode is characterized by a residue 
representing 80% of total system but multifragmcntation 
configurations arc already possible. For instance, config- 
urations characterized by two intermediate size fragments 
{Zmax ~ 30 and Z maX 2 « 20), though five times less 
probable than the most probable fragmentation mode, 
are nevertheless frequent enough to induce a quite flat 
Y(Z max ). The diversity of fragmentation modes trans- 
lated in broad Z max and Z max <x distributions persists at 
6 MeV/nucleon, but there is no more possible to identify 
a close competition among different fragmentation pat- 
terns. This means that there are no more distinct man- 
ners of filling up the available volume, whose co-existcncc 
leads to large fluctuation of the Coulomb energy. 

We remind at this point that fragmentation patterns 
are nevertheless very sensitive to break-up fragments def- 
inition or modelisation of the break-up stage itself. If, 
for instance, one sticks to the non-interacting break-up 
fragments scenario but considers that, in agreement with 
Thomas-Fermi calculations, excited nuclei at freeze-out 
arc diluted, the fragment charge distribution will be set- 
tled by the competition between the reduced free volume 
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FIG. 3: (Color online) Charge distributions (top) and distri- 
butions of average kinetic energy as a function of fragment 
charge (middle, bottom) corresponding to the multifragmen- 
tation of the (130,60) nucleus with V = 6Vb and E ex — 4 
MeV/nucleon. The two upper panels depict break-up stage 
results, while after-burner data are illustrated in the bottom 
panel. Solid circles stand for all fragments while open circles, 
open squares and open triangles stand for the largest, second 
largest and, respectively, third largest fragment in each event. 
The solid line on the top panel corresponds to the break-up 
fragment charge distribution obtained under the assumption 
that break-up fragments interact not only through repulsive 
hard-core and Coulomb potentials, but also via proximity po- 
tentials. 



and the augmented thermal kinetic energy. The same 
qualitative situation is reached if, not the fragments den- 
sity but, their internal excitation is modified. If, for ex- 
ample, one adopts for the nuclear level density an expres- 
sion which leads to lower fragment internal excitation, in 
view of Eq. ([5]) Kth will increase, favoring an increased 
reaction products multiplicity. This last quantity, in its 
turn, by making possible a more uniform population of 
the freeze-out volume characterized by a larger Vc, will 
tend to diminish Kth- 

Much dramatic modifications are expected if one con- 
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FIG. 4: (Color online) The same as in Fig. [3] for the mul- 
tifragmentation of the (130,60) nucleus with V = 6Vo and 
E ex = 6 MeV/nucleon. 



siders that break-up fragments interact not only via re- 
pulsive Coulomb, but also via attractive nuclear prox- 
imity potentials. This conceptually different approach is 
mainly justified by the fact that for break-up volumes 
of the order of few Vq the distances between fragment 
surfaces may be inferior than ~ l fm. This situation 
has been discussed recurrently by Das, De, Samaddar, 

' [3 111 together 
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Satpathy, Bonasera and collab. 
with break-up fragment subsequent recombination and 
shown to lead to an increased productivity of light and 
heavy fragments at the cost of the intermediate ones. As 
recombination, which occurs if two fragments approach 
each other during the Coulomb propagation, acts in the 
sense of washing-out the statistical properties of break- 
up fragment formation, here we shall restrict ourselves 
to comment exclusively on the consequences of modify- 
ing fragments energetics. 

It is relatively easy to anticipate from Eq. {5) that, by 
considering an extra attractive potential, one will get an 
increase of thermal kinetic energy and reaction products 



multiplicity. The confirmation is given by the solid curves 
on the top panels of Figs. [3] and [4] obtained in the case 
in which the nuclear interactions are implemented as in 
Ref. 
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In both situations one may notice a dramatic 
enhancement of the light cluster multiplicity and the to- 
tal supression of fragments with Z > 10. These steep 
Y(Z) distributions and the evolution of their slopes with 
source excitation may be reconciled with experimental 
data if and only if one assumes that final fragment forma- 
tion is dominated by post-break-up dynamics (including 
collective flow) and multiparticle correlation [12j . If this 
were the case, the freeze-out would occur much later that 
the break-up. The complete modclisation of this process 
is nevertheless a challenging task which goes beyond the 
goal of the present paper. 

In addition to charge distributions, fragment average 
kinetic energy distributions represent robust and directly 
accessible experimental information and make up a key 
ingredient in the standard procedure of identifying the 
statistically equilibrated source by confrontation with 
predictions of statistical models [8|, [9| . The middle and 
bottom panels of Figs. [3] and [4] depict the average ki- 
netic energy distributions of primary and, respectively, 
cold fragments corresponding to the same source (130,60) 
with V = 6Vo and E ex = 4 and 6 MeV/nucleon. The in- 
clusive distributions are plotted with full circles, while 
distributions corresponding to the largest, second largest 
and third largest fragment are plotted with open circles, 
squares and triangles. Collective radial flow is set to zero 
to keep fragments statistical properties unaffected. 

Fragment average kinetic energy distributions are qual- 
itatively similar for the two source excitations. As one 
may notice, the maximum value of 49 (42) MeV reached 
by the primary (asymptotic) < K(Z) > | 4 M eV/micieon 
distribution exceeds by 25% (30 %) the maximum value 
obtained by the corresponding < K(Z) > | 6 MeV/nucleon 
distribution. This result, in apparent contrast with what 
one would expect given the increase with 57% of the total 
Coulomb energy over the considered energy domain, may 
be understood taking into account the much stronger in- 
crease in the total number of reaction products [18[ . 

A common and interesting feature is present in the 
charge domain where the < K(Z) > distributions reach 
their maximum. Thus, the break-up and asymptotic av- 
erage kinetic energies of the largest fragment are sys- 
tematically smaller than the average kinetic energies of 
the second largest fragment which are, in their turn, 
smaller than the ones corresponding to the third largest 
fragment. This result has been already pointed out 
by the Indra collaboration in the case of Xe+Sn at 32 
MeV/nucleon and Gd+U at 36 MeV/nucleon reactions 
[lit an( l shown to diminish with the source excitation 
energy, in perfect agreement with the present results. 
Taking into account that fragment kinetic energies are 
to a large extent dictated by Coulomb, it becomes ob- 
vious that analyzing them one may get information on 
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FIG. 5: (Color online) Radial probability distributions of dif- 
ferent size (Z=l, 5, 10, 20, 30, 40 and 50) primary fragments 
at break-up. The statistically equilibrated source (130, 60) is 
characterized by an excitation energy of 4 MeV/nucleon and 
a freeze- out volume V = 6Vo- 



the most probable fragment position at break-up. Hav- 
ing the same dependence as Coulomb on fragment mass 
and distance from the freeze-out volume center, radial 
collective flow, if present, would enhance this shift. If 
this reasoning is correct, it means that larger a fragment 
is, closer to the freeze-out volume center it is produced. 

The answer to this issue is offered by Fig. [5] where ra- 
dial probability distributions of different size fragments 
corresponding to the break-up stage of the (130,60) 
source with V = 6Vq and E ex = 4 MeV/nucleon arc 
plotted. As a first general remark, one may say that the 
probability to create a fragment inside the freeze-out vol- 
ume, whatever its size, is highly un-uniform and strongly 
diminishes in the core region. Moreover, heavy fragments 
are localized preferentially towards the inner parts, while 
relatively light nuclei may be created over wider regions. 
This means that lighter is a fragment, stronger will be the 
Coulomb repulsion the charged core will exercise over it 
and, consequently, higher its final kinetic energy. This 
explains the observed systematic shift between the max- 
imum values of kinetic energy corresponding to the three 
largest fragments. The systematic reduction of the vol- 
ume accessible to a fragment as its mass increases is the 
consequence of the employed non-overlapping condition 
between a fragment and the wall of the container which 
mimics the freeze-out volume. In the case of the heaviest 
fragments (Z=50), this geometric condition is responsi- 
ble for fragment concentration in a region which repre- 
sents only 15% of the total freeze-out volume. We remind 
that the classification of multifragmentation events with 
respect to fragments' spatial arrangement and its influ- 
ence on fragment-fragment correlation functions was dis- 
cussed for the first time in Ref. [2(| in the framework 
of MMMC model Q , where the authors have identified 



FIG. 6: (Color online) Total charge radial distributions cor- 
responding to the (130,60) multifragmentating nucleus with 
V — 6Vb and different excitation energies (4, 6, 8 and 10 
MeV/nucleon). 

'sun'- and 'soup'-like events. 

As the source excitation energy increases and fragmen- 
tation becomes more advanced, a more uniform popula- 
tion of the freeze-out volume is expected, such that the 
largest fragments kinetic energy shifts become negligi- 
ble. The evolution of the total charge radial distribution 
with source excitation is illustrated in Fig. [6] for the 
same multifragmcnting nucleus, (130, 60) with V = 6Vq. 
Indeed, at 8 MeV/nucleon the matter in the inner re- 
gions of the freeze-out volume is 10 times denser that 
the one produced at 6 MeV/nucleon, but the overall dis- 
tribution remains strongly outwards peaked, giving the 
source a bubble-like structure. Bubble-like structures of 
the nuclear matter at break-up have been obtained also 
in the framework of stochastic mean-field approaches plj 
which explain fragmentation on behalf of growing volume 
and surface instabilities encountered during the expan- 
sion phase of the excited system, as recently reported in 
Refs. [22|, . This agreement between results of sta- 
tistical models with cluster degrees of freedom and dy- 
namical models with nucleonic degrees of freedom is far 
from being trivial taking into account the conceptually 
different scenarios the two categories of models advance 
for explaining multifragmentation and the almost com- 
plementary treatment of the physical process. 



V. CONCLUSIONS 

To conclude, using a microcanonical multifragmenta- 
tion model with cluster degrees of freedom we have ana- 
lyzed the break-up fragmentation patterns of a medium 
size equilibrated source who follows different paths 
through the excitation energy - freeze-out volume space. 
The constraints imposed on the freeze-out volume are 
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found to not affect significantly the magnitude of differ- 
ent energy fluctuations. Moreover, kinetic and configu- 
rational energy fluctuations are insensitive to the system 
phase properties as far as the considered fragment par- 
titions arc similar. Over the whole domain of excitation 
energy, spatial matter distribution at break-up is highly 
un-uniform, its outward peaked shape giving the source 
a bubble-like structure. The most probable localization 
of nuclear fragments at break-up depends on fragment 
mass and, because of Coulomb acceleration, it is possible 
to infer it from the experimentally accessible fragment 
average kinetic energy distributions, especially at inter- 
mediate values of source excitation. Thus, heavy frag- 
ments are found to be produced in the inner regions of 
the freeze-out volume, while the lighter ones are produced 
in a larger region of the freeze-out volume. Considering 
that break-up fragments interact not only through repul- 
sive hard-core and Coulomb potentials but also via prox- 
imity potentials, one obtains dramatic modifications of 
the break-up fragmentation patterns which suggest that 
final fragment formation is strongly influenced by post- 
break-up dynamics and multiparticlc correlations. 
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